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Abstract 



A variety of phenomena in nuclear and high energy physics seemingly do 
not satisfy the basic hypothesis for possible stationary states to be of the 
type covered by Boltzmann-Gibbs (BG) statistical mechanics. More specifi- 
cally, the system appears to relax, along time, on macroscopic states which 
violate the ergodic assumption. Some of these phenomena appear to follow, 
instead, the prescriptions of nonextensive statistical mechanics. In the same 
manner that the BG formalism is based on the entropy Sbg = ~~ k J2iPi l n Pi> 
the nonextensive one is based on the form S q = k(l — J2iPi)/( a ~ 1) (with 
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Si = Sbg)- Typically, the systems following the rules derived from the former 
exhibit an exponential relaxation with time toward a stationary state char- 
acterized by an exponential dependence on the energy (thermal equilibrium), 
whereas those following the rules derived from the latter are characterized by 
(asymptotic) power-laws (both the typical time dependences, and the energy 
distribution at the stationary state). A brief review of this theory is given here, 
as well as of some of its applications, such as electron-positron annihilation 
producing hadronic jets, collisions involving heavy nuclei, the solar neutrino 
problem, anomalous diffusion of a quark in a quark-gluon plasma, and flux of 
cosmic rays on Earth. In addition to these points, very recent developments 
generalizing nonextensive statistical mechanics itself are mentioned. 



I. INTRODUCTION 

The foundation of statistical mechanics comes from mechanics (classical, quantum, rel- 
ativistic, or any other elementary dynamical theory). Consistently, in our opinion, the 
expression of entropy to be adopted for formulating statistical mechanics (and ultimately 
thermodynamics) depends on the particular type of occupancy of phase space (or Hilbert 
space, or Fock space, or analogous space) that the microscopic dynamics of the system under 
study collectively favors. In other words, it appears to be nowadays strong evidence that 
statistical mechanics is larger than Boltzmann-Gibbs (BG) statistical mechanics, that the 
concept of physical entropy can be larger than 

w 

Sbg = -k^pilnpi (1) 

i 

(hence Sbg — k In W for equal probabilities), pretty much as geometry is known today to be 
larger than Euclidean geometry, since Euclid's celebrated parallel postulate can be properly 
generalized in mathematically and physically very interesting manners. 

Let us remind the words of A. Einstein [1] expressing his understanding of Eq. (1): 
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Usually W is put equal to the number of complexions [...]. In order to calculate W , one needs 
a complete (molecular-mechanical) theory of the system under consideration. Therefore it 
is dubious whether the Boltzmann principle has any meaning without a complete molecular- 
mechanical theory or some other theory which describes the elementary processes. S = 
jjlogW + const, seems without content, from a phenomenological point of view, without 
giving in addition such an Elementartheorie. 

This standpoint is, in our opinion, quite similar to the position that Riemann adopted 
when he began his study of "the concepts which lie at the base of Geometry" . Along this 
line, it is our understanding that the entropy expressed in Eq. (1) is, on physical grounds, no 
more irreducibly universal than Euclidean geometry with regard to all possible geometries, 
which, by the way, also include the fractal one, which inspired the theory addressed in this 
paper. Nonextensive statistical mechanics [2], to which this brief review is dedicated, is 
based on the following expression 

1 _ yw q 

S « = k q Z\ ( q e ^ ^ = S bg) ■ (2) 

It is well known that for microscopic dynamics which relax on an ergodic occupation of 
phase space, the adequate entropic form to be used is that of Eq. (1). Such assumption 
is ubiquitously satisfied, and constitutes the physical basis for the great success of BG 
thermostatistics since over a century. We believe that a variety of more complex occupations 
of phase space may be handled with more complex entropies. In particular, it seems that Eq. 
(2), associated with an index q which is dictated by the microscopic dynamics (and which 
generically differs from unity), is adequate for a vast class of stationary states ubiquitously 
found in Nature. In recent papers, E.G.D. Cohen [3] and M. Baranger [4] also have addressed 
this question. 

A significant amount of systems, e.g., turbulent fluids ( [5,6] and references therein), 
electron-positron annihilation [7,8], collisions of heavy nuclei [9-11], solar neutrinos [12,13], 
quark-gluon plasma [14], cosmic rays [15], self-gravitating systems [17], peculiar velocities 
of galaxy clusters [18], cosmology [19], chemical reactions [20], economics [21-23], motion of 
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Hydra viridissima [24], theory of anomalous kinetics [25], classical chaos [26], quantum chaos 
[27], quantum entanglement [28], anomalous diffusion [29], long-range- interacting many- 
body classical Hamiltonian systems ( [30] and references therein), internet dynamics [31], 
and others, are known nowadays which in no trivial way accommodate within BG statistical 
mechanical concepts. Systems like these have been handled with the functions and concepts 
which naturally emerge within nonextensive statistical mechanics [2,32,33]. 

We may think of q as a biasing parameter: q < 1 privileges rare events, while q > 1 privi- 
leges common events. Indeed, p < 1 raised to a power q < 1 yields a value larger than p, and 
the relative increase p q /p = p q ~ l is a decreasing function of p, i.e., values of p closer to (rare 
events) are benefited. Correspondingly, for q > 1, values of p closer to 1 (common events) 
are privileged. Therefore, the BG theory (i.e., q — 1) is the unbiased statistics. A concrete 
consequence of this is that the BG formalism yields exponential equilibrium distributions 
(and time behavior of typical relaxation functions), whereas nonextensive statistics yields 
(asymptotic) power-law distributions (and relaxation functions). Since the BG exponential 
is recovered as a limiting case, we are talking of a generalization, not an alternative. 

To obtain the probability distribution associated with the relevant stationary state (ther- 
mal equilibrium or metaequilibrium) of our system we must optimize the entropic form (2) 
under the following constraints [2,32]: 

£ft = 1 > (3) 

i 

and 

^iPi Ei = JJ U\ 

V rfl q ' 1 ' 

t—ti fi 

where {Ei} is the set of eigenvalues of the Hamiltonian (with specific boundary conditions), 
and U q is a fixed and finite number. This optimization yields 

Pi — y 5 K°) 

where 



z q = £[i - (i - q )(3 q ( Ej - cgf/a-s) , (6) 

3 

and 

/5 being the optimization Lagrange parameter associated with the generalized internal energy 
U g . Equation (5) can be rewritten as 

Pi oc [1 - (1 - 9)Wd-«> ee ef E * , (8) 

where j3' is a renormalized inverse "temperature" , and the q- exponential function is defined 
as e x = [1 + (1 — g)^ 1 ^ 1- ^ = 1/[1 — (g— l)^] 1 /^" 1 ) (with = e 21 ). This function replaces, in a 
vast number of relations and phenomena, the usual BG factor. In particular, the ubiquitous 
Gaussian distribution oc e~ ax2 becomes generalized into the distribution oc e~ aqx2 = 1/[1 + 
(q - l)a q x 2 Yl^ (fat-tailed if q > 1). 

II. GENERALIZING NONEXTENSIVE STATISTICAL MECHANICS 

Nonextensive statistical mechanics generalizes the BG theory. It presumably addresses 
(multi) fractal-like occupancy of phase space at the associated stationary states (e.g., metae- 
quilibrium), instead of the usual, homogeneous, occupancy which satisfies ergodicity. Is 
there any fundamental reason for stopping here? We do not think so. In fact, roads point- 
ing towards generalizations of (or alternatives for) nonextensive statistical mechanics are 
already open in the literature. Let us mention here two of them (already exhibiting some 
successes), namely (i) crossovers between g-statistics and g'-statistics ( [15] and references 
therein), and (ii) the recently introduced Beck-Cohen superstatistics [34]. Both of them 
address the energy distributions corresponding to the stationary states, and are perfectly 
compatible, as we shall show. More precisely, the first type can be thought as a particular 
case of the second type. However, statistical mechanics is much more than a distribution 
correctly corresponding to the stationary state. Indeed, if any given entropy S({pi}) is op- 
timized by a specific distribution p iy all entropic forms which are monotonic functions of 
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S will be optimized by the same distribution. Nevertheless, only a very restricted class of 
such entropic forms can be considered as serious candidates for constructing a full thermo- 
statistical theory, eventually connected with thermodynamics. In particular, one expects 
the correct entropy to be concave and stable. Such is the case [35] of S q as well as of the 
generalized entropy recently proposed [36,37] for the just mentioned superstatistics [34]. We 
briefly address these questions in this Section. 

Let us first consider the following differential equation: 

d f x =ay (y(0) = 1). (9) 

The solution is given by 

y = e ax . (10) 

We can use this result in at least three manners which are of interest in statistical mechanics: 

(i) We may refer to the sensitivity to the initial conditions, and consider x = t, where t is 
time, y = £ = liniA;r(o)->o Ax(t) / Ax(0) , where Ax(t) denotes the discrepancy of two initial 
conditions in a one-dimensional map (or, for higher- dimensional systems, the analogous 
situation for the axis along which the maximum dynamical instability occurs), and a = 
Ai 7^ 0, where Ai is the Lyapunov exponent. In this case Eq. (10) reads in the familiar form: 

m=e Xlt . (11) 

(ii) We may refer to the relaxation towards the stationary state (thermal equilibrium), 
and consider x = t, y = [0(t) — O(oo)]/[O(0) — 0(oo)], where O is the average value of a 
generic observable, and a = — 1/r < 0, where r is a relaxation time. In this case Eq. (10) 
reads in the typical form: 

9^. - - e -t/r (m 

O(0)-O(oo)~ 6 ■ [ ] 

(iii) We may refer to the energy distribution at thermal equilibrium of a Hamiltonian 
system, and consider x = Ei, where Ei is the energy of the i-th microscopic state, y = 
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Zp(Ei), where p is the energy probability and Z the partition function, and —a = (5 is the 
inverse temperature. In this case Eq. (10) reads in the familiar BG form: 

p-PEi 

P{E l ) = — (Z = 5>-^)- (13) 

Z 3 

This distribution is of course the one that optimizes the entropy Sbg under the standard 
constraints for the canonical ensemble. 

Let us next generalize Eq. (9) as follows: 

d f x = *y q (2/(0) = 1). (14) 

The solution is given by 

y = e^ ! = [l + (l-q)ax] 1 ^ , (15) 

e q being from now on referred to as the g-exponential function (ef = e x ). The three above 
possible physical interpretations of such solution now become 

(i) For the sensitivity to the initial conditions, 

e(*)=ej«* = [l + (l-9)\*] 1/(1 " 9) , (16) 

where \ q ^ is the generalized Lyapunov coefficient (see [26]), and, at the edge of chaos, 
X q > and q < 1. 

(ii) For the relaxation, 

°(*) ~ = e -Vr = 1 (17 ) 

C(0)-C(oo) q [1 + (q - l)r,*]V(«-i) ' 1 ; 

where r 9 > is a generalized relaxation time, and typically q > 1 [38]. 

(hi) For the energy distribution, we get the form which emerges in nonextensive statistical 
mechanics, namely [2,32] 

- ^ = (z « s V E,) • (i8) 

where usually, but not necessarily, f3' q > and g > 1. This distribution is the one that 
optimizes the entropy S q under appropriate constraints for the canonical ensemble. 



Let us next unify Eqs. (9) and (14) as follows: 



dy_ 

dx 



aw + a q y q (y(0) = 1) . 



(19) 



This is a simple particular case of Bernoulli equation, and its solution is given by 



V = 



e (i- q )a 1 x_ 1 \f 1 + 1 
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(20) 



This solution reproduces Eq. (10) if a q = 0, and Eq. (15) if a\ = 0. It corresponds to a 
crossover from a q ^ 1 behavior at small values of x, to a q — 1 behavior at large values of 
x. The crossover occurs at x c ~ l/[(<? — [38]. 



Let us now briefly review five recent applications of the ideas associated with nonextensive 
statistical mechanics to phenomena in nuclear and high energy physics, namely electron- 
positron annihilation [7,8], collisions of heavy nuclei [9-11], the solar neutrino deficit [12,13], 
quark anomalous diffusion [14], and the flux of cosmic rays [15]. 
Electron-positron annihilation: 

In high energy collisions of an electron with a positron, annihilation occurs and, immedi- 
ately after, typically two or three hadronic jets are produced. The probability distribution 
of their transverse momenta is non-Boltzmannian, more strongly so with increasing energy 
of collision. This phenomenon has defied theoreticians since several decades, particularly 
since Hagedorn [16] quantitatively analyzed such data in the frame of BG statistical me- 
chanics. A phenomenological theory has been recently proposed by Bediaga et al [7] , which 
beautifully fits the data. The fitting parameters are two, namely the temperature T and the 
entropic index q. For each energy collision E c a set of (T, q) is determined. It numerically 
comes out that q depends on the energy (like q(oo) — q(E c ) oc E' 1 / 2 for increasingly large 
E c , and q(E c ~ 0) ~ 1), but T does not! This invariance of T with respect to the colli- 
sion energy constitutes the central hypothesis of the physical scenario advanced long ago by 



III. APPLICATIONS 
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Hagedorn. This scenario is now confirmed. The ingredients for a microscopic model within 
this approach have also been proposed [8]. 
Heavy nuclei collisions: 

A variety of high-energy collisions have been discussed in terms of the present nonexten- 
sive formalism. Examples include proton-proton, central Pb-Pb and other nuclear collisions 
[9,10]. Along related lines, entropic inequalities applied to pion-nucleon experimental phase 
shifts have provided strong evidence of nonextensive quantum statistics [11]. 
Solar neutrino problem: 

The solar plasma is believed to produce large amounts of neutrinos through a variety 
of mechanisms (e.g., the proton-proton chain). The calculation done using the so called 
Solar Standard Model (SSM) results in a neutrino flux over the Earth, which is roughly 
the double of what is measured. This is sometimes referred to as the neutrino problem or 
the neutrino enigma. There is by no means proof that this neutrino flux defect is due to 
a single cause. It has recently been verified that neutrino oscillations do seem to exist ( 
[12] and references therein), which would account for part of the deficit. But it is not at 
all clear that it would account for the entire discrepancy. Quarati and collaborators [13] 
argue that part of it - even, perhaps, an appreciable part of it - could be due to the fact 
that BG statistical mechanics is used within the SSM. The solar plasma involves turbulence, 
long-range interactions, possibly long-range memory processes, all of them phenomena that 
could easily defy the applicability of the BG formalism. Then they show [13] in great detail 
how the modification of the "tail" of the energy distribution could considerably modify the 
neutrino flux to be expected. Consequently, small departures from q — 1 (e.g., \q — 1| of 
the order of 0.1) would be enough to produce as much as 50% difference in the flux. This 
is due to the fact that most of the neutrino flux occurs at what is called the Gamow peak. 
This peak occurs at energies quite above the temperature, i.e., at energies in the tail of the 
distribution. 
Quark diffusion: 

The anomalous diffusion of a charm quark in a quark-gluon plasma has been analyzed by 
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Walton and Rafelski [14] through both nonextensive statistical mechanical arguments and 
quantum chromodynamics. The results coincide, as they should, only for q = 1.114. 

Cosmic rays: 

The flux $ of cosmic rays arriving on Earth is a quantity whose measured range is 
among the widest experimentally known (33 decades in fact). This distribution refers to a 
range of energies E which also is impressive (13 decades). This distribution is very far from 
exponential: See Figs. 1 and 2. This basically indicates that no BG thermal equilibrium is 
achieved, but some other (either stationary, or relatively slow varying) state, characterized 
in fact by a power law. If the distribution is analyzed with more detail, one verifies that 
two, and not one, power-law regimes are involved, separated by what is called the "knee" 
(slightly below 10 16 eV). At very high energies, the power- law seems to be interrupted 
by what is called the "ankle" (close to 10 19 eV) and perhaps a cutoff. The momenta 
Mi = ((E-(E)) 1 ) = [fo cutoff dE §(E)(E — (E)) l ]/[fQ CUtoff dE ${E)] {I = 1,2,3) as functions 
of the cutoff energy E cuto ff (assumed to be abrupt for simplicity) are calculated as well: See 
Figs. 3, 4 and 5. At high cutoff energies, (E) saturates at 2.48944... x 10 9 eV [15], a value 
which is over ten times larger than the Hagedorn temperature (close to 1.8 x 10 s eV [8]). In 
the same limit, we obtain for the specific-heat-like quantity M 2 ~ (E 2 ) ~ 6.29 x 10 21 (eV) 2 . 
Finally, M 3 ~ (E 3 ) diverges with increasing E cuto ff- This is of course due to the fact that, 
in the high energy limit, $ oc 1/E~~ 2 ~ 1/E 3A ; consequently the third moment integrand 
vanishes like 1/E 0A , which is not integrable at infinity. 

One may guess that, along such wide ranges (of both fluxes and energies), a variety of 
complex intra- and inter-galactic phenomena are involved, related to both the sources of 
the cosmic rays as well as the media they cross before arriving on Earth. However, from a 
phenomenological viewpoint, the overall results amount to something quite simple. Indeed, 
by solving a simple differential equation, a quite remarkable agreement is obtained [15]. This 
differential equation is the following one: 

|| = -b' P i - b P j . (2i) 
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This differential equation has remarkable particular cases. The most famous one is (q 1 , q) = 
(1,2), since it precisely corresponds to the differential equation which enabled Planck, in 
his October 1900 paper, to (essentially) guess the black-body radiation distribution, thus 
opening (together with his December 1900 paper) the road to quantum mechanics. The 
more general case q' = 1 and arbitrary q is a simple particular instance of the Bernoulli 
equation, and, as such, has a simple explicit solution (Eq. (20) with a± = —b' and a q = —b). 
This solution has proved its efficiency in a variety of problems, including in generalizing the 
Zipf-Mandelbrot law for quantitative linguistics (for a review, see Montemurro's article in 
the Gell-Mann-Tsallis volume [33]). Finally, the generic case q > q' > 1 also has an explicit 
solution (though not particularly simple, but in terms of two hypergeometric functions; 
see [38]) and produces, taking also into account the ultra-relativistic ideal gas density of 
states, the above mentioned quite good agreement with the observed fluxes. Indeed, if we 
assume < b' « b and q' < q, the distribution makes a neat crossover from a power- 
law characterized by q at low energies to a power-law characterized by q' at high energies, 
which is exactly what the cosmic rays exhibit to a quite good approximation. Let us finally 
mention that the first possible microscopic interpretation of our phenomenological theory 
has just been suggested [39]. 

For possible effects of a slightly nonextensive black-body radiation on cosmic rays see 
[40]. Finally, other aspects related to cosmic rays have been shown to exhibit fingerprints 
of nonextensivity [41]. 

IV. CONCLUSIONS 

In nuclear and high energy physics, there is a considerable amount of anomalous phe- 
nomena that benefit from a thermostatistical treatment which exceeds the usual capabilities 
of Boltzmann-Gibbs statistical mechanics. This fact is due to the relevance of long-range 
forces, as well as to a variety of dynamical nonlinear dynamical aspects, possibly leading 
to nonmarkovian processes, i.e., long-term microscopic memory. Some of these phenomena 
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appear to be tractable within nonextensive statistical mechanics, and we have illustrated 
this with a few typical examples. For the particular case of cosmic rays, we have indi- 
cated their average energy (E) ~ 2.48944... x 10 9 eV, and the specific- heat-like quantity 
(E 2 ) — (E) 2 ~ 6.29 x 10 21 {eV) 2 , with the hope that they might be usefully compared to 
related astrophysical quantities, either already available in the literature, or to be studied. 
Along the same vein we have also presented the dependence of such momenta on a possibly 
existing high-energy cutoff. 

In addition to this, we have sketched the possible generalization of nonextensive statistical 
mechanics on the basis of a recently introduced entropic form [36], whose stationary state 
is the Beck-Cohen superstatistics [34]. The metaequilibrium distribution associated with 
a crossover between g-statistics and ^'-statistics can be seen as a particular case of this 
generalized nonextensive statistical mechanics. 

It is worthy to mention at this point that the present attempts for further generalization 
of BG statistical mechanics are to be understood on physical grounds, and by no means as 
informational quantities that can be freely introduced to deal with specific tasks, and which 
can in principle depend on as many free (or fitting) parameters as one wishes. Examples 
of such informational quantities are the Renyi entropy (depending on one parameter and 
being usefully applied in the multifractal characterization of chaos), Sharma-Mittal entropy 
(which contains both Renyi entropy and S g as particular cases), and very many others that 
are available in the literature. The precise criteria for an entropic form to be considered a 
possible physical entropy are yet not fully understood, although it is already clear that it must 
have a microscopic dynamical foundation. It seems however reasonable to exclude, at this 
stage, those forms which, in contrast with Sbg and S q , (i) are not concave (or convex), since 
this would seriously damage the capability for thermodynamical stability and for satisfactory 
thermalization of different systems put in thermodynamical contact, or (ii) are not stable, 
since this would imply that the associated quantity would not be robust under experimental 
measurements. These crucial points, and several others (probably equally important, such as 
the finite entropy production per unit time), have been disregarded by Luzzi et al [42] in their 
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recent criticism of nonextensive statistical mechanics. Indeed, Renyi entropy, Sharma-Mittal 
entropy (that Luzzi et al mention without any justification at the same epistemological level 
as S g ) are neither concave nor stable for arbitrary values of their parameters. These and 
other information measures (most of them not concave and/or not stable) have been freely 
introduced, along various decades, as optimizing tools for specific tasks. They can certainly 
be useful for various purposes, which do not necessarily include the specific one we are 
addressing here: a thermodynamically meaningful generalization of the Boltzmann-Gibbs 
physical entropy. 
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FIG. 1. Flux of cosmic rays as function of their energy. See [15] for details. 
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FIG. 2. Same as in Fig. 1, but the ordinate is now multiplied by (Energy) 3 . 

3 




T 



1 




< 




1 





1x10 



10 



5x10 
E 



n 



lxl0 12 -| 



cutoff 



io 15 

'cutoff f eV ] 

FIG. 3. (E) function of the cutoff energy. 
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FIG. 4. {E 2 ) (black dashed curves) and M2 = (E 2 ) — (E) 2 (red continuous curves) as functions 
of the cutoff energy. 
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FIG. 5. {E 3 ) (black dashed curves) and M 3 = (E 3 ) -3{E)(E) 2 + 2{E) 3 (red continuous curves) 
as functions of the cutoff energy. At vanishing E cuto fj, M3 vanishes from below, i.e., with slightly 
negative values. 
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